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A study of the radiative neutrino counting reaction e + e~ — ► vvj at LEP1 and LEP2 
energies is presented. An approximate expression for the spectrum of the observed 
^ . photon is derived within the framework of the pt-dependent structure function ap- 

proach. This is compared with an exact expression and found in agreement within 
the foreseen experimental accuracy. This model describing single-photon radiation 
can be applied to the more general case of initial-state single-photon emission accom- 
panying invisible final-state events. Higher-order QED corrections due to undetected 
initial-state radiation are also included. The implementation in a Monte Carlo event 
generator is briefly described. 
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1 Introduction 



The radiative neutrino-counting reaction and, more in general, the single- 
photon reactions, provide a very useful tool for measuring the number of light 
neutrinos and eventually detecting new-physics signals at LEP 0. 

The lowest-order spectrum for single-photon production within the Stan- 
dard Model is known [0] . For energies around the Z° peak this exact calculation 
agrees within about 1% with two different approximations ||, |J . 

Concerning higher-order QED corrections, two different calculations and 
the corresponding codes do exist ||, [|. Both are based on the QED structure 
function approach ||. In the first one M, QED corrections are implemented on 
the tree-level spectrum in the so-called Point Interaction Approximation (PIA) 
and the corresponding program is a semi-analytical code. In the second one ||, 
0(a) plus higher-order QED corrections are included for the Z°-exchange only 
and the contribution of the W^-exchange diagrams is added at the tree level. 
This formulation is implemented in a Monte Carlo generator of unweighted 
events. However, as is known, the treatment of the QED radiative corrections 
around the Z° peak can become inadequate as the energy increases. A proper 
treatment of the hard photonic radiation has to be provided. 

The purpose of this work is twofold: 

i) to improve the approximations existing in the literature for the study of the 
reaction at LEP1 energies, in order to get reliable results for radiative neutrino 
counting as well as for single-photon events at LEP2; to provide a new result 
on the treatment of the hard-photon spectrum; 

ii) to establish a general strategy for developing a single-photon library, inclu- 
sive of all the standard and non-standard signals. This last point is motivated 
by the relevance of the experimental studies of single-photon radiation accom- 
panying invisible events at LEP for the search of possible new-physics effects. 

The outline of the paper is as follows. In Sect. 2 the exact and approxi- 
mate calculations for the spectrum of the radiative neutrino counting reaction 
at lowest order in the Standard Model are reviewed. A quantitative study 
of the approximations at LEP2 versus LEP1 energies will be performed and 
an approximate expression for the photonic spectrum will be presented. In 
Sect. 3 it will be shown how to include in the formulation higher-order initial- 
state QED corrections, and our results will be compared with existing ones. 
Some single-photon distributions of experimental interest will be shown and 
commented in Sect. 4. Finally, conclusions and possible perspectives on the 
subject are given in Sect. 5. 
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2 Lowest-Order Calculations 



The production of neutrino-antineutrino pairs and of a single photon is de- 
scribed within the Standard Model of the electroweak interactions by s- and 
t-channel diagrams involving Z° and W ± exchange Jl|. W contributions occur 
only for the production of electron neutrinos and they can be distinguished in 
diagrams with single-VT exchange and diagrams with the exchange of two W's, 
the latter involving the non-Abelian trilinear coupling WW'y. On the basis of 
the above diagrams, the differential cross section in the photon energy and 
angle cos $ 7 has been calculated exactly and its expression contains the Z° 
and W squared amplitudes and the W-Z° interference. What is important to 
stress about this calculation is that it is performed keeping a finite value for 
the mass of the W boson and, therefore, it applies without restrictions on the 
energy range. 

However, for the energies around the Z° peak, two approximations are 
known in the literature and they both reproduce, with an accuracy of the 
order of 1%, radiative neutrino counting experiments at LEP1. The two ap- 
proximations are: 

• The Point Interaction (or Contact) Approximation [R| (PIA). It consists 
in neglecting the diagram with exchange of two W's and taking the limit 
Myj/ — > oo in the calculation. 

• The Point Interaction Approximation with Angular Radiator Q (PIAH). 
It uses the QED ^-dependent structure functions, i.e. the method con- 
taining also the transverse degrees of freedom of the emitted radiation [0. 
It consists in dressing the cross section for producing a vv pair (with- 
out photons) with an angle- dependent radiator. The photon spectrum is 
given in this approximation by the following factorized expression: 

= a£ IA ((l - x 7 )a)H^(x^a). (1) 

Lli Jb ry LL ry 

a? IA is the cross section for e + e —^vv as evaluated within the PIA ap- 
proximation at the energy scale reduced by photon emission; (X~, Cry] S 
is the angle- dependent radiator which describes the probability of radi- 
ating a photon with a given energy fraction x 1 = E^/Ef, at the angle 
$ 7 (c 7 = cos$ 7 ) . Its expression, derived from the 0(a) ^-dependent 
structure function based on the electron — > electron + photon splitting 
function, reads as M 
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where P(l — x) is the Altarelli-Parisi splitting function and pt e -,+ rep- 
resent the electron and positron transverse momenta, respectively. In 
terms of the photon angle c 7 , one obtains || 



a 1 + (1 — x 7 ) 2 
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This expression contains the leading infrared and collinear singularities 
represented by the l/x 7 term and by the angular factor in (fj), respec- 
tively. This specific expression can, however, be properly modified to 
include less dominant (next-to-leading) contributions, as will be shown 
in the following (see [f|). 

Figure 1 shows the integrated radiative neutrino counting cross section 
with the minimum photon tagging angle = 20° and the minimum photon 
energy E™ m = 1 GeV. In this plot, the solid line is the result obtained by 
integrating the exact photon spectrum of Ref . [Q] , the dotted and dashed lines 
correspond to the PIA |§ and PIAH Q approximations, respectively. As 
can be seen, the two approximations are in good agreement with the exact 
calculation around the Z° peak, but they start to differ as the energy goes 
above about 120-130 GeV and the difference is well visible at higher LEP2 
energies. Figure 2 is a blow-up of Fig. 1 at LEP1 (Fig. 2a) and LEP2 (Fig. 2b) 
energies respectively. The two approximations agree with the exact calculation 
within ~ 1% at LEP1 (where they essentially provide the same integrated cross 
section with a hard photon tail above the peak due to the emission of a Z°- 
return photon of energy E 1 ~ (1 — M\j s) y/s/2), whereas at LEP2 they are 
in net disagreement with the full spectrum. In particular, they provide an 
increase in cross section with energy, as a consequence of the absence of finite 
VU-mass effects needed at LEP2 energies. Furthermore, the larger deviation 
of the PIAH approximation is due, in particular, to the specific expression 
employed for the angular radiator as in eq. @. 

However, hard photons, strongly suppressed around the Z° peak, are im- 
portant at LEP2 and they have to be included. Within the framework of the 
Pt-dependent structure function approach, by virtue of its generality, a simple 
approximate expression can be achieved, which constitutes one of the results 
of the present paper. This is of some importance for experiments at LEP2 be- 
cause, if a proper expression for handling the photon spectrum is found, it can 
be used as a basic tool to develop a library of single-photon events, including 
standard and non-standard (in particular SUSY) processes. The approximate 
expression can be implemented with the two following steps: 

• to replace the cross section e + e~ — > vv in PIA approximation (<Jq IA ) with 



3 



the corresponding exact cross section in the Standard Model (<jQ Xact ), i.e. 
for e + e~ — > (Z, W) — > uu, thus including the necessary Vi^-mass effects (]; 

to adopt, instead of the form in eq. (0), the more accurate expression 
for the angular radiator in such a way that also a part of next-to-leading 
contributions are included. This properly modified radiator has been 
already derived M and reads: 



rr{a) a 1 
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The latter expression differs from the radiator of eq. (FJ) by less than 1% around 
the Z° peak, but it is more appropriate to describe the emission of hard photons 
(of energy fraction x 1 ~ 1) at large angles. 

Therefore the photonic spectrum of the reaction e + e~ — > vv^ can be writ- 
ten in factorized form as follows: 



d-P _ pcxat i ,, ; ,.\ u"' ] 

dec ryd/C? 



((1 - x 7 )s)ifS(x 7 ,c 7 ;s). (4) 



Let us comment about the physical meaning of the above equation. Given 
the specific bare kernel e + e~ — ► uu, as well as for any other kernel correspond- 
ing to any given process of the type e + e~ invisible objects, dressing it with 
the angular radiator Hg^ amounts to attaching a photon line on the external 
charged legs; in this way, the "universal" , factorized form of the photonic ra- 
diation will be included, comprehensive of some of the next-to-leading effects. 
Therefore, in this approximation the W diagram with the photon emitted by 
the internal propagator is not reproduced. 

However, by virtue of the generality of the method, if the cross section for 
an e + e~ scattering yielding an invisible final state is known, an approximation 
of the cross section for the corresponding radiative process can be obtained in 
a straightforward way. 

Figure 3 shows a comparison of the exact spectrum with the approximate 
results obtained using the two different angle- dependent radiators of (0) and 
(HD, respectively. The solid line corresponds to the exact spectrum; the dotted 
and dashed lines are the result derived from convoluting the exact bare cross 
section with the angular radiators of eq. @ and eq. respectively. 

As can be clearly seen, whereas the approximation based on the radiator of 
eq. (0) differs by several per cent from the exact result, using the "modified" 
radiator of eq. (§) improves the agreement at the level of a few per cent. To 
be more precise, for 150 < ^/s < 175 the relative difference between the exact 

4 The exact invisible cross section was computed by means of SCHOONSCHIP H. 
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results and the approximation of eq. ([J) is of the order of 1-2% and becomes of 
the order of 3-4% only at about 200 GeV. These deviations, depending on the 
centre-of-mass energy, can be seen as a measure of the non-leading effects not 
included in the angular radiator approximation. However, for the experimental 
precision foreseen at LEP2, a theoretical error of the order of a few per cent 
is more than adequate. 



3 Higher-Order QED Corrections 



To match the experimental precision expected at LEP2, the corrections due to 
multiphoton soft emission and to the radiation of hard photons emitted in the 
very forward direction (and hence lost in the beam pipe) by the colliding beams 
must be included. The large QED corrections introduced by the undetected 
initial-state radiation can be successfully described within the framework of 
the QED structure function approach. This method, as recalled in the Intro- 
duction, is indeed employed in both the formulations developed for radiative 
neutrino counting at LEP1 M, |5|. However, both the approaches have their 
drawback and therefore need to be updated for experiments at LEP2. A sim- 
ple solution, which represents an improvement of both the existing approaches 
and has an overall accuracy of about 1%, can be achieved by convoluting 
the exact spectrum of Ref. || with electron and positron structure functions. 
The QED-corrected cross section can then be written, in analogy with QCD 
factorization, as a convolution of the following form f9j: 



o (s 



[ dx\ dx2 dE 1 dc y D(x\, s)D(x2, s) — , (5) 
J dE^dCy 



where da/dE^dcj is the exact spectrum of e + e~ — > i/u^y |J, the photon vari- 
ables being referred to the centre-of-mass frame after initial-state radiation, 
and D(x,s) is the electron (positron) structure function. Its expression, as 
obtained by solving the Lipatov-Altarelli-Parisi evolution equation in the non- 
singlet approximation, is given by M: 



+ —I3 2 
32 M 



r (i + 
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D(x,s) = IV n |(l-#-i-^(l + x) (6) 
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where L = In (s/m 2 ) is the collinear logarithm, 7^ is the Euler constant 
and r(z) is the gamma function. The first exponentiated Gribov-Lipatov 
term describes multiphoton soft emission, the second and third ones hard 
bremsstrahlung in the collinear approximation. Further non-leading QED cor- 
rections are included as a .fT-factor, computed from the exact 0(a) result 



of Ref. ||10fl . The relevant weak corrections are taken into account, for the 
Z°-exchange contribution, in the form of Improved Born Approximation, ac- 
cording to the general recipe of Ref. ]ll] . 



The master formula (gD is implemented in a Monte Carlo event genera- 
tor |TJ| Q. Figure 4a shows the QED-corrected cross section as a function of 
the centre-of-mass energy at LEP2. The dash-dotted line represents the exact 
cross section at lowest order; the solid line is the result corresponding to eq. (^j) 
(i.e. in the case of convolution of the full spectrum), and the dotted line, re- 
ported for the sake of comparison, shows the results obtained by simulating the 
approach of Ref. H], namely correcting the Z° contribution only, and adding to 
this result the TV-exchange diagrams at tree level. Two considerations are in 
order here. First, the QED-corrected cross section is higher than the Born one 
as a consequence of the Z° radiative return: the effect is to enhance the Born 
result by a factor of about 1.3. Secondly, the convolution of the full spectrum 
is in good agreement (within 1%) with the approach of Ref. f| because the 
QED-corrected cross section is largely dominated by the Z° radiative return, 
and the tree-level contribution of W diagrams (and W-Z interference) is al- 
most flat over the full energy range spanning from LEP1 to LEP2 (Fig. 4b). Of 
course, a more general (beyond the Standard Model) single-photon spectrum, 
including initial-state QED corrections, can be simply obtained. This can be 
achieved by substituting, in eq. (5), the exact proper generalized spectrum (or, 
if sufficient, its approximation according to eq. (4)) for da/dE^dc^ within the 
Standard Model. 

4 Simulation of Photon Distributions 

In this section, some technical details about the Monte Carlo event generator 
developed on the basis of the formulation presented in the previous section 
will be sketched; a complete and exhaustive account of the code will be given 



elsewhere |fL3 |. A sample of realistic experimental distributions obtained by 
analysing the single-photon events generated by the Monte Carlo code will also 
be shown and briefly discussed. 



5 For a first attempt to develop a neutrino counting Monte Carlo based on p t -dependent 
structure functions, see Ref. [fl2|. 
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Since the integrand is a superposition of different peaking behaviours, the 
standard importance-sampling technique, within the multichannel approach, 
was employed to develop an efficient event generator. The QED-corrected 
cross section is written as the sum of four different channels with the following 
features: (1) the radiative variables x 12 are generated according to the soft 
exponentiated part of the structure function, and the photon energy is then 
generated according to the dependence of the Z° propagator on E 7 ] (2) x 12 
as above and E 1 generated according to the infrared-like behaviour l/E^; 
(3) E 1 flat, x\ as the soft part of the structure function and x 2 according 
to the dependence of the Z° propagator on the reduced centre-of-mass energy 
s = X1X2S; (4) the same as (3) with the interchange of x\ <-> x 2 . Moreover, to 
take care of the Lorentz boost caused by the emission of unbalanced electron 
and positron radiation, the photon phase-space variables are first generated 
in the centre-of-mass frame and then boosted to the laboratory frame by the 
boost given by the Lorentz factor (3 L = {x\ — x 2 )j (x 1 + x 2 ) valid for collinear 
kinematics. 

Figures 5, 6 and 7 show the photon energy distribution, the photon angu- 
lar distribution and the pj / E b distribution at a fixed LEP1 (E b = 48 GeV) 
and LEP2 (E b = 88 GeV) beam energies, respectively, assuming an experi- 
mental apparatus with E 1 — 1 GeV and = 20°. The solid histograms 
correspond to the lowest-order approximation, the dashed ones include the ef- 
fect of initial-state radiation. In order to consistently compare the tree-level 
and QED-corrected distributions, the numbers of events integrated in the two 
cases are proportional to the corresponding integrated cross sections. Con- 
cerning the photon energy distribution (Fig. 5), two peaks are well visible 
both at LEP1 (Fig. 5a) and LEP2 (Fig. 5b) energies: the higher one is located 
at the energy value of about (I — M|/s) v /i/2, the lower one is due to l/E^ 
peaking behaviour. As can be seen, the main modifications introduced by 
initial-state radiation are to reduce the higher peak and to enhance the lower 
one, as expected by typical convolution effects, thus modifying the two-bump 
asymmetrical shape of the Born distribution into a more symmetrical one. The 
photon angular distribution (Fig. 6) has the same symmetrical shape (peaked 
in the forward direction) at LEP1 (Fig. 6a) and LEP2 (Fig. 6b) energies, the 
only difference being that at LEP2 the QED-corrected distribution is higher 
than the Born one because of the Z° radiative return. At last, comparing 
the Pt/Eh distributions at LEPI (Fig. 7a) and LEP2 (Fig. 7b) energies the 
emission of mostly soft photons at LEPI and of harder ones at larger angles 
at LEP2 can be clearly recognized. 



7 



5 Conclusions 



In this paper a study of the neutrino counting reaction e + e~ — > uuj at LEP2 
energies has been presented. Within the ^-dependent structure function ap- 
proach, an approximate expression for the photon energy and angle spectrum 
has been derived. The accuracy of the approximation, as estimated by com- 
parison with the exact result, is of the order of 1-2% near the V^-pair threshold 
and grows larger (up to about 3-4%) only at about 200 GeV. By virtue of its 
generality, the approximation can be used as a first attack strategy to dress 
standard as well as non-standard (e.g. SUSY) processes giving rise to an in- 
visible final state in order to get an approximate formula for the spectrum 
of the corresponding radiative single-photon events. Therefore, on the basis 
of these considerations, a library of single-photon processes should be feasible 
for experiments at LEP2. A further investigation on the tree- level generalized 
spectra would be of course required for processes whose form is less simple 
than the one shown in eq. (4). 

Concerning higher-order QED corrections, a formulation based on structure 
function convolution of the exact photon spectrum is presented, which can be 
estimated to be accurate at the 1% level and agrees with the approximate 
approach of Ref. |J within about 1%. This is implemented in a Monte Carlo 



event generator [13], which also can be easily generalized to describe single- 



photon radiative processes accompanying invisible events at LEP. 
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Figure Captions 



Figure 1. The lowest-order integrated cross section of e + e~ — > z/z/7 for E 1 = 
1 GeV and = 20° as a function of the centre-of-mass energy. Solid line: 
exact spectrum; dotted line: PIA; dashed line: PIA with angular radiator. 

Figure 2. Blow-up of Fig. 1 at LEP1 (Fig. 2a) and LEP2 (Fig. 2b) energies. 

Figure 3. The lowest-order integrated cross section at LEP2. Solid line: exact 
spectrum; dotted line: "soft" angular radiator of eq. (2); dashed line: "modi- 
fied" angular radiator of eq. (3). 

Figure 4. The QED-corrected integrated cross section at LEP2 (Fig. 4a). 
Solid line: convolution of the exact spectrum; dotted line: convolution of the 
Z° contribution with W at tree level; dash-dotted line: exact Born spectrum. 
The W and W-Z° interference contributions as a function of the centre-of- 
mass energy (Fig. 4b). 

Figure 5. The photon energy distributions at LEP1 (Fig. 5a) and LEP2 
(Fig. 5b) energies. Solid histogram: Born; dashed histogram: QED-corrected. 

Figure 6. The same as Fig. 5 for the photon angular distribution. 

Figure 7. The same as Fig. 5 for the pi/Eb distribution. 
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